Heatmap any DEG

library(pheatmap)
data = read.delim('~/DEGs_Infected_vs_Control_Combined.csv')

anycommondeg = data[data$Number_UP != 0 | data$Number_DOWN != 0,3:7]

anydeg = data[data$Number_UP >= 2 | data$Number_DOWN >= 2, 3:7]

rowd = c("correlation","euclidean", "maximum", "manhattan", "canberra", "binary", "minkowski")
cold = c("correlation","euclidean", "maximum", "manhattan", "canberra", "binary", "minkowski")
hclust  = c("ward.D", "ward.D2", "single", "complete", "average", "mcquitty" , "median"  , "centroid")

for (i in rowd){
    for (k in hclust) {
      pheatmap(as.matrix(anydeg[,2:5]), 
               col = topo.colors(100), 
               clustering_distance_rows = i,
               clustering_distance_cols = i,
               clustering_method = k,
               main = paste0(i,k),
               labels_col = c('Adults','Children','Mice - Feet','Mice - LN'))
    }
}